Granular flow down a rough inclined plane: transition between thin and thick piles 
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The rheology of granular particles in an inclined plane geometry is studied using molecular 
dynamics simulations. The flow-no-flow boundary is determined for piles of varying heights over a 
range of inclination angles 8. Three angles determine the phase diagram: 8 r , the angle of repose, 
is the angle at which a flowing system comes to rest; 8 m , the maximum angle of stability, is the 
inclination required to induce flow in a static system; and 8 ma x is the maximum angle for which 
stable, steady state flow is observed. In the stable flow region 8 r < 8 < 8 ma x, three flow regimes can 
be distinguished that depend on how close 8 is to 8 r : i) 8 >> 8 r : Bagnold rheology, characterized 
by a mean particle velocity v x in the direction of flow that scales as v x oc h 3 ^ 2 , for a pile of height 
h, ii) 8 > 8 r : the slow flow regime, characterized by a linear velocity profile with depth, and iii) 
8 ~ 8 r : avalanche flow characterized by a slow underlying creep motion combined with occasional 
free surface events and large energy fluctuations. We also probe the physics of the initiation and 
cessation of flow. The results are compared to several recent experimental studies on chute flows 
and suggest that differences between measured velocity profiles in these experiments may simply be 
a consequence of how far the system is from jamming. 
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I. INTRODUCTION 

The flow properties of granular materials puzzle engineers and physicists alike. Despite numerous experimental 
studies of heap flows, avalanches, and chute flows, there still does not exist a general consensus on the nature of 
typical granular flow. Early systematic experiments in inclined plane geometries were concerned with flow rates of 
granular material down chute , usually with smooth bases. Large regions of plug flow were observed, i.e. the whole 
bulk of the material slid down the plane at a velocity determined by the input mass rate and inclination angle of the 
chute. Kinetic theories for the inclined plane problem can explain some aspects of this flow behaviour for nearly elastic 
particles on smooth chutesBi. Other studies have determined structural features in the flow in two dimensional (2D) 
experiments! and simulations!'!. Still, there has not been the same success in predicting flows down rough inclines 
for inelastic, frictional particles. 

Detailed measuring techniques have been developed recently to capture various flow properties of granular mate- 
rial down rough inclined planesll~0 and in heap flow geomctries@~0. Accompanying these more thorough experi- 
ments, further theoretical treatments have also emergecflllUli, together with full three dimensional (3D) computer 
simulations&il. Despite this growth in literature, there appears to be a disturbing feature of these studies: reports 
on quantities such as velocity profiles differ quite dramatically. While some reports are agreement with the pioneering 
results of Bagnold§10, others are not. 

In 1954, Bagnoldil proposed a momentum transfer picture to describe his experimental results on incrtial granular 
flow. This simple 2D model describes the collisional process between particle layers during flow. The resulting relations 



predict a shear stress a that is proportional to the square of the strain-rate, 



er oc 7 , 



(1) 



where the strain-rate is defined as the derivative of the velocity in the direction of flow (x) with respect to depth (z), 
7 = dv x (z)/dz. Typically, in local, Newtonian, hydrodynamic treatments of granular dynamicsS, the relationship 
a = 777 is used, combined with the kinetic term r\ oc T 1 / 2 , where 77 is a viscosity and T is the granular temperature (to 



be discussed in section III B). Since the strain-rate is also proportional to the square root of the granular temperature 
7 oc T 1 / 2 , Eq. ^ is then recovered. Bagnold's argument, applied to the case of bulk granular flow, is predicated on a 
constant density depth-profile with no-slip condition at the base. This leads to a velocity depth-profile of the form, 



v x {z) oc \h 3/2 ~{h-z) 



3/2 



(2) 



for a pile of total height h, where z is measured from the bottom of the pile. Because Bagnold scaling comes naturally 
from dimensional analysis in the absence of intrinsic time scales, one might expect that most studies would observe 
this kind of rheology. However, results depend strongly on the experimental procedure; profiles of this form have only 
been reported in a few instances!. 

One method to study gravity-driven flows is to construct a wedge-shaped static pile that evolves to produce an 
inclined surface close to the angle of repose of the materiaJll 0, Fig. |l|(a). We term these heap flows. Here, the 
inclination of the surface is determined by material properties. To study flow, material is supplied via a hopper-feeder 
mechanism at the high end of the pile and flow is induced at the free surface. In this case, although the flow rate may 
be varied by changing the discharge rate of the hopper, the system typically flows only in a thin surface layer, indicated 
by the shaded region in Fig. |l](a), which is of order 10 particles diameters in depth. In one such experiment, Lemieux 
and DurianEl observed avalanche behavior for low particle flux, characterized by sporadic downstream movement 
at the free surface. On increasing the particle flux, a velocity profile with depth was observed to be approximately 
exponential in the thin surface layer. More recently, utilizing a similar experimental geometry, Komatsu et alEl found 
that the 'stationary' bulk material beneath the surface flow actually undergoes very slow, 'creeping' motion with a 
well-defined exponential velocity profile at longer times. 





FIG. 1. Schematics of experiments used to study gravity-driven, granular flows. The (a) heap flow geometry used by Lemieux 
and DurianEl and Komatsu et alEl, and (b) the inclined plane geometry employed by Pouliquen for chute flow studies!! Our 
simulations correspond closely to the Pouliquen experiments. 



Another method uses inclined planes to study flows and avalanches, typically with rough beds designed to reduce 
slip at the base, thereby avoiding plug or unstable flowEBilS. We term these chute flows. Roughness is introduced 
onto the chute by gluing particles to the chute base!0@, or by attaching a frictional material to the base, such as 
the velvet cloth used by Daerr and DouadyEH. Chute flow studies of two and three dimensional systems with rough 
bases report depth-averaged velocities with both linear00 and non-linear dependenc<J30 on height. The geometry 
for these experiments is shown in Fig. |l](b). The material is contained in a reservoir at one end of a plane inclined at 
a variable angle 9 with respect to the horizontal. After the gate to the reservoir is opened, material flows down the 
incline. Steady flow occurs only at high enough angle; while flows come to rest quickly on the rough bed for small 
inclination anglesEI Increasing the gate opening leads to thicker flows down the incline: there are no dead zones 
along the incline and all the material flows down the chute. This allows one to independently control the inclination 
angle and material flux of the pile, in contrast to the heap experiments. 

Using a range of inclination angles and flow depths, PouliquenEI measured the depth-averaged flow velocity of 
material for piles of different thickness. The velocity in the direction of flow scaled with height of the pile as v x cx ft 3 / 2 , 
which has a different form from the heap flow results reported by Lemieux and Durianlll and Azanza et al.0. However, 
Ancey's comprehensive experimental investigation0 demonstrated that during the course of one experiment, a range 
of flow behaviors is obtained that depends, among other parameters, on the flow rate of the material. Ancey's results 
suggest that the different flow profiles observed are not necessarily in contradiction, but may be a consequence of the 
conditions under which the measurements are made. 

In our earlier studies0il, we simulated the rhcology of dense, cohesionless, frictional spheres in inclined plane 
geometries. Employing periodic boundary conditions in the flow and vorticity directions, our simulations can be 
considered equivalent to the experimental set up in Fig. 0(b), for an infinitely long chute far away from the material 
reservoir and side walls, with a continuous supply of material. Our previous study was concerned with the flow of 
thick piles on rough bases over a range of tilt angles where we found Bagnold scaling in stable steady flow, for a range 
of parameters (e.g. particle friction and inelasticity, and inclination angles). 

To investigate the discrepancies between the various experimental results previously mentioned, we have extended 
our previous computational studies of granular chute flows to cover a wide range of pile thickness, from very thin 
to very thick piles, both near and far from the angle of repose 9 r . We show here that even thick piles experience a 
crossover in behavior from Bagnold rheology to linear velocity profiles, and also avalanche-like flow near the angle 
of repose. We therefore refine our definition of stable, steady state flow: we report on three distinct steady state 
flows observed in our simulations which match the experimentally observed flows reported by Lemieux and DurianEl 
(avalanches), Ance (linear velocity profiles), and Pouliqucri (Bagnold scaling). Our main conclusion is that the 
observed behavior is strongly dependent on how close the system is to jamming or how close 9 is to 9 r . 



In the next section we briefly describe the simulation procedure. In section III we present the results of our 
simulations, in particular defining the phase space of our study and describing transitions between the flow properties 
of thick and thin piles both near and far from the angle of repose. We also present results for our analyses on the 
initiation and cessation of flow in section In section |v| we summarize the main conclusions of this work. 



II. SIMULATION PROCEDURE 



Here wc briefly introduce the simulation model; a thorough description of the technique is available elsewhereEJEi 
We model N cohesionless, frictional, and inelastic spheres of diameter d and mass m flowing on a fixed, rough bed of 
area A with a free top surface. Although most of our results are for monodisperse systems, we have verified that our 
results hold for systems with small polydispersityS 

We employed a 3D contact force model that captures the major features of granular particles: interactions between 
spheres are modeled as Hertzian contacts with static friction. The implementation of the contact forces, both the 
normal forces and the shear (friction) tangential forces, is essentially a reduced version of that employed by Walton 
and Braun0, developed earlier by Cundall and StrackH. Normal and tangential forces acting on particle i due to 
contact with particle j are given by 



F ™; 3 = V S v/ d {knSijiiij - m e ff7„v ny ) , (3) 
F Uj = ySij/d (-k t u tij - m e g7 t v t(3 .) , (4) 

where Sij is the particle-particle overlap, = r^/r^- = (r^ — i*j)/|rj — r*j- j , is the unit normal along the line of centers, 
for particle i situated at Yi and j at Yj. v n> t are relative velocities in the normal and tangential directions respectively, 
and Uf y is the elastic tangential displacement that is set to zero at the initiation of a contact. m c ff is the effective 
mass and k Ut t and j ni t are elastic and viscoelastic constants respectively. Thus in a gravitational field g, the force on 
particle i arises from Newton's second law for the translational and rotational degrees of freedom, 

Ff t = mg + 2 F »«+ F *«; ( 5 ) 

j 

and the force on particle j is determined from Newton's third law. 

As shown in Fig. |^, we define flow to be parallel to the ^-direction, z the height normal to the base, and y, the 
vorticity axis, is normal to the xz-plane. Values for the friction coefficient fi, and the elastic and viscoelastic coefficients 
(which are related to the moduli of the material) are chosen to be characteristic of real materials. For all simulations 
in this study we set /i = 0.50, k n = 2 • 10 5 mg/d, k t — |fc„, 7„ = 50^/gJd, and jt — 7^/20. The time-step St = 10~ 4 t, 
where r = \J d/g. For a millimeter sized particle r ss 0.01s. All simulations were run for at least 1 • lO 7 ^ before any 
quantities were measured. 

Figure shows snapshots of flowing systems where the area of the chute base A = Lx x Ly = 20rf x lOd = 200d 2 , 
and we vary the total number of particles N. The fixed chute bed is made from a random close packing of particles 
with the same diameter as those in the bulk, mimicking gluing particles onto a chute base as in Ref.i, for example. 



In section III B we kept A = 200eP fixed. In section [V, we used a base with dimensions A — 200d x lOd — 2000cZ 2 . 
We have shownEl that results for thick piles are robust with respect to the dimensions of the chute base. The height 
of a static pile h — N/Ap depends on the density of the static packing p 0l or equivalently the volume fraction 
(f> = Po tt/6 (r> 0.600) which is defined in the hard sphere limit (where the subscripts denote the static values) . In the 
flowing state, the actual height of the pile increases with inclination angle due to dilation effects of the flow. However, 



the dependence of h, and therefore p, on 8 is weak. We find that the computed bulk volume fractions vary over a 
small range, 0.55 < <f> < 0.60, for steady state flows for the particular parameters chosen in this paperEl. Thus for the 
simulations h Q = ttN/6A(/> : for N = 5000, h 22d. 




(c) (d) 

FIG. 2. Configuration snapshots for systems of various sizes with a chute base of fixed area A = 200d 2 , that is inclined at an 
angle 6 = 24° to the horizontal, where (a) N = 8000, (b) N = 6000, (c) N = 4000, and (d) N = 2000. The fixed rough base is 
denoted by the black particles, the flowing particles are grey, and the frame is a guide to the eye fixed at the height h ~ 39. 6d 
of the N — 8000 system to provide a comparison between the different flow heights. 



III. RESULTS 

In our previous studies0jH, we roughly determined the range of parameters where steady state flow existed for 
chute flows, focusing mainly on steady state flow for thick piles away from 9 r . We found that 9 r is approximately 
independent of height for piles of sufficient depth. In this work we perform a more thorough analysis of the entire region 
of the h — 8 phase space. We determine several different flow regimes within the stable, steady state region which are 
qualitatively identical to several experimental observations in inclined plane geometries, on avalanche dynamicsEl'0 



and rheolog Although the exact location of the boundaries of the phase picture shown in Fig. |^ are dependent 

on specific particle parameters, such as the friction and damping coefficients, the features we report here are general 
to systems that exhibit steady state flow. 



A. Phase Behavior 



The overall phase behaviour is determined by the height of the pile h and the angle of inclination 9, as characterised 
in 2D by Pouliquen and Renault!!. A prominent feature of that study is the existence of a flow line, later observed 
in 3D experiments by Pouliquen! and also in studies of avalanchesEl. The flow line defines the angle of repose as 
a function of height of material flowing 9 r (h). One can equivalently refer to this in terms of the conjugate variable 
/istop(#)& 9 r (h) is the angle at which a pile of thickness h will cease to flow, and conversely, h s t op {9) gives the 
maximum thickness of material that will cease to flow at an inclination angle 9. Various theoretical treatments have 
also successfully predicted the existence of the flow lineB!!. Experimental!!'!! and numerical!! studies have also 
identified distinct sub-regions within the stable flow regime of which more will be discussed below. 

In order to quantify this boundary with precision for our system, we performed a number of simulations to generate 
the h — 9 phase diagram shown in Fig. || Initially, a static packing with a volume fraction <f> w 0.60, at 9 = O°0, is 
tilted until flow is observed. The maximum angle of stability 9 m is defined as the angle at which flow is initiated and 
largely depends on construction history of the initial static packing, such as the roughness of the supporting bed and 
the packing fraction of the initial configuration!!. The angle of repose 9 r is defined as the inclination angle at which 
the system subsequently jams and flow ceases upon reducing 9 back down. 
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FIG. 3. The h — phase diagram for granular chute flow, where h is the height of the pile. The filled squares are the 
simulation data for 6 r (h) and we indicate the flow line as a fit to Eq. |^. The open squares are 8 m (h). Solid circles are 6 max (h), 
the angle at which we observe the onset of unstable flow. Route A corresponds to a path taken by keeping the angle fixed and 
reducing the system size N, whereas Route B has N fixed and decreasing inclination. These routes are discussed in section 



[HE. 



Our initial packing was generated by taking a flowing system and gradually reducing the inclination down to 0° and 
then allowing the system to relax until it became completely static. Experimental studies on avalanches have shown 



that 9 m can fluctuate over a range of 3° — 4° for different packings with the same height because of construction 
history effectsEl. 

We have also observed avalanche-like dynamics in our simulations. We characterise avalanches as distinct from 
regular or continuous flow by observing the time evolution of bulk quantities of the system such as bulk kinetic 
energy or stress. For continuous flow, the stress or bulk energy will fluctuate about a well-defined mean value, with 
fluctuations typically less than 10% of the mean. In the case of avalanches, the stress and energy fluctuate about a 
mean value with 0(1) fluctuations, indicative of sporadic flow that mainly occurs at the free surface. We also checked 
any possible system size dependence. To compare with our smaller systems, we studied a system with N — 32000 and 
A = 80d x 40d (i.e. approximately 10d high on a base with an area four times our standard system size). We again 
observed avalanche-likc behaviour for the big system, although we cannot rule out that in the limit of infinite system 
size, avalanches might not be observed. 

From our extensive study, we were able to determine the following general features of flow down inclined planes: 
there are three principal regions, corresponding to no flow, stable flow, and unstable flow. For a pile of given thickness 
h, and fixed microscopic interaction parameters, these three regions are separated by two angles: 9 r , the angle of 
repose, and max , the maximum angle for which we observe stable flow. 

No flow 9 < 9 r : the system does not flow. Starting from 9 = 0°, an inclination of 9 — 9 m > 9 r is required to 
initiate flow. As 9 m is approached, a few particles at the free surface rearrange and tumble short distances 
before coming to rest, but this does not lead to bulk motion or flowH. 9 r is the angle at which flow ceases on 
reducing the inclination from a flowing state - it is the lowest angle at which flow is observed and depends on 
material properties (friction and inelasticity )@, as well as the roughness of the supporting bed. 

Stable flow 9 r < 9 < 9 max : bulk flow is observed (after flow is initiated). The bulk averaged kinetic energy of the 
system is constant with time and the entire pile is seen to be in motion while no slip is observed at the bottom 
with these rough boundary conditions. The stable flow region can be divided into the following three regimes: 

Continuous flow I 9 >> 9 r : away from the flow line, the rheology is well approximated by Bagnold scaling, 
Eqs. ® and @). 

Continuous flow II 9 r < 9 < 9 m : close to the flow line, the rheology is not Bagnold- like, rather the velocity 
profile is better approximated by a linear relationship with depth v{z) oc z. The upper boundary 9 m is also 
seen in experimental studies of avalanches and differentiates between different flow mechanisms, depending 
on whether one is above or below m B 

Avalanching flow 9^9^: Very close to the angle of repose, but just above it, regular continuous flow is 
not observed, rather in this regime, steady state involves slow, underlying creep motion of the entire pile 
combined with sporadic free surface motion down the incline. In the vicinity of 9 r , extremely large energy 
fluctuations are observed0. 

Unstable flow 9 > 9 max : the development of a shear thinning layer at the bottom of the pile results in lift-off and 
unstable acceleration of the entire pile. 9 max becomes approximately constant for h > lOd. 

The flow descriptions above capture the principal observations of several different experiments of gravity-driven 
nowsi'SEIil. In particular, the transition from no flow to avalanche flow at 9 sa 9 r is similar to the observations 



tan#. r = tan 0i + [tan 62 — tan#i] exp 



(7) 



reported by Lemieux and DurianE3 for low material flux. The transition from avalanches to our continuous flow II 
regime appears to resemble the observations reported by Daerr and Douadj@ of avalanche studies for 9 r < 9 < 6 m . 
Similar observations have been reported for rotating cylinder experiments^. 

PouliquerJI, later verified by Daerr and Douadyil, characterized the flow line - the dividing line between no flow 
and flow. Using an empirical argument, Pouliquen suggested the following relation: 

h 

7 

where 9\ is the angle at which 9 r becomes independent of depth, 62 is the largest value of 9 r , and I is a characteristic 
length scale. In Fig. |§| we fit Eq. ^ to our simulation points and obtained the following fitting parameters; 6% = 
19.60° (19.55°), 6*2 = 28.50° (27.75°), and I « 8.5d, where the values in parentheses are the data points from the 
simulations. Our simulations not only predict the phase behavior of chute flows, but also capture the subtle features 
emerging in recent experiments. In the following sections we present various depth profiles of density, velocity, and 
stress to better quantify the transitions between thin and thick piles. 



B. Depth Profiles 



We computed depth profiles as a function of z for the packing fraction 0, velocity in the direction of flow v x , 
strain-rate 7, shear stress ct xz , and granular temperature T (see definition below), for a range of system sizes and 
inclination angles. Since our simulations employ periodic boundary conditions in the plane perpendicular to height, 
our simulations are not influenced by the effects of side walls and are therefore equivalent to measurements made 
along the centre line of a real chute flow, away from the edges. It has been recognised that profiles taken at the side 
walls are significantly different from those made along the centre lineQfl 

Because of the dependence of 9 r on h, we approached 9 r in two ways: fixing the inclination angle of the chute and 
observing flow in systems of different heights until flow ceased for a pile of height h stop (9) (route A). Alternatively, 
we held N fixed and reduced the inclination angle of a flowing state down to 9 r (h) (route B). 

Route A: Figures show various bulk profiles for route A, that is, systems at a fixed inclination angle of 9 = 24° 
but of varying sizes, N = 8000, 5000, 3000, 2000, 1000. The density profiles in Fig. | exhibit a constant density 
region throughout the bulk of the material which is a familiar observation for thicker systems!!. 




FIG. 4. Depth profiles of the packing fraction <j>, for systems of various size N, varying in heights from h m 40d to 5d for 
= 24°: route A in Fig. B. 



The flow velocity v x and strain-rate 7 profiles shown in Fig. |B| are for the same systems reported in Fig. For 
moderate to thick piles, h > 15d, the velocities and corresponding strain-rates scale with height. However, for h < 15d, 
there are significant differences. The velocity profiles, top panel in Fig. |^, for the thinner piles at first become more 
linear with depth, then ultimately have a curvature opposite to that of the thicker piles. This gradient transformation 
is more apparent from the strain-rate curves (bottom panel Fig. [5]): for thicker piles 7 decreases from the chute base 
toward the free surface of the flowing pile, whereas in thinner piles there is a range of system sizes where 7 ~ constant 
throughout the pile, leading to a linear velocity profile. For thinner systems still, 7 increases near the free surface. 




0.4 0.6 

FIG. 5. Depth profiles of the velocity in the direction of flow (top panel) and the corresponding strain-rate (bottom panel), 
each normalized by their maximum value, for various system sizes for 9 = 24° . Note the depth is normalized by the respective 
height for clarity. Legend as in Fig. W. 



In this work, we define the "granular temperature" T to be the averaged velocity fluctuations about the mean bulk 
flow velocity at that height, T = (v 2 ) — (v) 2 . The temperature in the flow direction is defined as T x — (v 2 ) — (v x ) 2 ■ 
Figure ^| shows the change in the profiles of T x from thick to thin piles, which is similar to the behavior of 7. From 
the simulations we find that for all sizes, 



T x oc 7 

although the data for thinner piles is considerably noisier. 

1.0 



(8) 




FIG. 6. Profile of the granular temperature T x /T™ ax , in the direction of flow as a function of z for various system sizes for 
= 24° . Legend as in Fig. §. 



Similarly, the rheology plot in Fig. ^| shows that thicker piles display Bagnold rheology: a xz oc j 2 , but there is a 
transition away from this as the piles become thinner. For thin piles, 5 < h < 15, the strain rate is approximately 
constant away from the free surface and the chute base. For the thinnest piles h « 5d, the strain rate decreases 
with shear stress. This may come about from competing length scales in the system, i.e. the thickness of the flow 
determines the separation between the top free surface and the chute base. Some theoretical treatments based on 
these considerations have been used to explain the differences between observed velocity profiles in flowing pilesB. 
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FIG. 7. Rheology data for various system sizes for 8 — 24°: j vs. a xz . Legend as in Fig. |4| 

Route B: Figures ||jll] show depth profiles for a system taken along route B in Fig. |3[ These are for a system of 
fixed size N = 6000, where we vary the inclination angle down toward 9 r . In Fig. ^| the density profiles for the system 
N = 6000 are shown for a range of inclination angles 9 = 25°, 23° , 22°, 21° , 20°. We see the constant density region, 
where the bulk density decreases with increasing angle due to dilatancy effects. However, there are no qualitative 
differences between the behaviour at different angles. 




FIG. 8. Depth profiles of the packing fraction 
inclination angles, route B in Fig. ^| 



for a system with a fixed number of particles, N = 6000, over a range of 



For the velocity and strain-rate profiles shown in Fig. ^, differences do emerge, and these differences mimic the 
behavior observed along route A. Away from the flow line, the velocity profiles scale as Eq. ||, whereas for systems 
at lower angles the velocity profiles change to linear. 7 changes likewise as does the granular temperature (shown in 



Fig. 0). 
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FIG. 9. Profiles for the velocity in the direction of flow v x (top panel) and the corresponding strain-rate 7 (bottom panel), 
each normalized by their respective maximum values, for a system with a fixed number of particles, N = 6000, over a range of 
inclination angles inclination angle. The height is normalized by the system size. Legend as in Fig. 0. 
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FIG. 10. Profile of the x-component of the granular temperature T x /T™ ax for a system with a fixed number of particles, 
N — 6000, over a range of inclination angles. Legend as in Fig. |§| 
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FIG. 11. Rheology data for a system with a fixed number of particles N 
normalized by (9 — 8 r ) vs. a xz . Legend as in Fig. H. 



6000, over a range of inclination angles: 7 



In Fig. [n] we plot the rheology curves for the different angles. We have rescaled 7 by (9 — 9 r ), for convenience. We 
find that for systems with (9 — 9 r ) > 1.0°, the rheology plots scale well. However, for the system that is closest to 9 r , 
the flow behavior is very different, though similar to the differences observed along Route A. 

Therefore the two routes A and B indicated in Fig. || offer complementary paths to similar states. Because of the 
very nature of the dependence of 9 r on h, the flow properties of systems at different angles and different heights show 
similarities provided they are taken equivalently close to 9 r . 



C. Velocity Scaling 

Pouliqueni showed that flows of thick and thin piles approximately obey a simple scaling relation. Non- 
dimensionalizing v, the velocity measured along the centre line at the top of the flowing pile by \fgh, where g is 
the acceleration due to gravity and h is the measured height of the pile, one obtains the Froude number Fr = v/y/gh. 
By plotting the Froude number against h/h stop , where h stop (9) is the height of a pile at which a system flowing at 
an angle 9 would stop, Pouliquen was able to collapse his full data set of systems with different particle size and 
properties. Pouliquen proposed the simple scaling relation, 

v 9 n n-stop 

where the coefficient f3 was fit to the experimental data, resulting in an apparently universal value, (3 exp = 0.136. 

The quantity h s t op encodes details of the particle roughness, the base roughness and all other intrinsic properties of a 
particular system. Accordingly, one need only measure h stop of the system in question to determine the flow properties. 
In a similar vein, in Fig. [l^ we plot the flow velocity, taken from the top of the flowing pile, non-dimensionalised by 
\J~gh for each height h against h/h s t op - The data agrees well with a linear relation, giving a coefficient p s i m — 0.147, 
which is in good agreement with the experimental result determined by Pouliqueni, (3 exp = 0.136. 
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FIG. 12. Froude number v/y/gR as a function of h/h s t op for different inclinations. A measure of the slope gives (3 S i m = 0.147, 
compared with Pouliquen's experimental determination of f3 exp — 0.136. 



The scaling relation Eq. g although simple, appears to be quite robust. Pouliquenc3 and also Andreotti et al. H. 
have taken Eqs. ^ and |?] to predict phase behaviour and flow properties. Their resulting analyses of granular fronts 
captures many of the main features of granular flows. However, as yet, there is no underlying explanation for these 
seemingly simple relationships. 



IV. INITIATION AND CESSATION OF FLOW 



We have also studied the initiation and cessation of flow, something difficult to probe experimentally. For these 
studies we used large systems, A = 200rf x lOd = 2000<i 2 and h a ~ 50d, providing better resolution for these processes. 
In addition, any size effects due to thin samples do not complicate our observations. 

We have previously studied initiation of flow from a crystalline static state in 2D@. After such a system is tilted, a 
disordering front propagates upwards from the bottom of the pile, melting the ordered structure, until the entire pile 
is disordered and flows in its steady state form. Although the starting state in 3D need not be ordered, we find that 
flow initiation in 3-D also undergoes a similar bottom-up process. 

Initiation was studied for a disordered, static packing of our larger system that was tilted from 8 — 0° to 24°. 
See Fig. [l3| At the instant of flow initiation, the system slips at the bottom boundary and begins to flow as a plug 
as indicated in Fig. |l3|(a). Particles closest to the base undergo collisional motion with the fixed bed leading to the 
creation of a granular temperature front, as shown in Fig. [T3|(b). This thermalization front propagates up the pile, 
disrupting the plug region. After the original initiation slippage, the slip region disappears and the velocity profile 
appears approximately linear with height, with a crossover to the plug region. After the plug phase melts, the system 
dilates slightly, leading to a velocity profile of the Bagnold steady state form. 




FIG. 13. Initiation of flow: time evolution of the depth profiles of, (a) the velocity in the direction of flow and (b) the 
granular temperature in the direction of flow, at the initiation of flow for a static state (h « 50d), initially at 9 — 0°, that is 
tilted to 24°. Initial failure occurs near the base, the pile then "heats up" from the bottom of the pile, a thermalization front 
propagates up the pile with time, disrupting the plug flow above it. The lowest curve starts at t = 025t after the system is 
first tilted, and we plot curves for t — 0.50r, It, 2t, 3t, and 4r respectively. 



Cessation was observed for a system that was initially in steady state flow at 9 — 24°, as the inclination angle 
was reduced to 6 — 18°, which is below 9 r . See Fig. |l4|. Although it may seem intuitive to expect an exact reversal of 
the initiation process, this is not observed. After 9 is reduced to below the angle of repose, particles at the bottom of 
the pile, which have the smallest velocities, come to rest first becoming trapped by the rough base. The frozen layer 
then serves to trap the slowly moving particles above it. This trapping front propagates up the pile towards the free 
surface, as indicated in Fig. |l4|(a). A consequence of this is that the granular temperature cools from the chute base 
as the whole system gradually comes to rest. We note the emergence of a thin saltating region at the free surface 
during the cooling process. 




FIG. 14. Cessation of flow: time evolution of the depth profiles of, (a) the velocity in the direction of flow and (b) the 
granular temperature in the direction of flow, for the cessation of flow of a system with h ~ 50d. Initially in steady state flow 
with = 24° , the inclination angle is then reduced to 18° < 9 r leading to the gradual stopping of flow into a jammed state. 
Particles moving slowest at the base stop flowing first, leading to a freezing front that propagates up the pile. The top curve is 
taken at a time lOOr after the angle is first reduced to below the angle of repose, and subsequent curves are separated by lOr 
later. 

During this trapping process the pile exhibits two distinct behaviors separated by the trapping front. Below the 
trapping front, the particles are frozen into a pack, whereas above the trapping front, particles continue to flow. As 
the trapping front moves up the pile, the depth of the flowing section continuously shrinks with time. As shown in 
Fig. |l4|(a), this leads to a bulk velocity profile that consequently transforms from Bagnold-like flow to linear with 
depth until flow ceases. Essentially, flow cessation is a transient version of route B in Fig. || for flow at one angle over 
a range of flow heights. 

Another feature of flow transitions is in understanding the stress response of the system upon initiation and cessation 
of flow. Classical treatments of granular medial! demand that failure will only occur if the ratio of the shear stress 
to the normal stress equals the tangent of the angle of repose, 

<J XZ = cr 22 tan6» r . (10) 

However, Eq. Il^ needs to be modified to take into account the history dependence of the starting point of flow for 



initiation at 9 m , and the general stopping point of flow for cessation at 9 r . In Fig. [15| the stress ratio o xz j a zz is plotted 
as a function of inclination angle 9 for initiation of flow (circles) and cessation of flow (squares). We see that a xz j 'a zz 
is linear with tan9. Still, to observe flow from initiation, the stress ratio must exceed tan# m , or for flow to cease its 
value must be reduced to below tan# r . In other words, the shear stress can have a non-zero value and remain static as 
long as the stress ratio is smaller than the tangent of 9 r or 9 m . The particular stress state of the static system is very 
much history dependent and comes about from the imposed conditions at the instant the system became static0. 
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FIG. 15. Ratio of the shear stress to the normal stress a xz /a zz as a function of the tangent of the angle of inclination 9 for 
initiation of flow (squares) and cessation of flow (circles). The open symbols correspond to static states and the filled symbols 
are flowing. The lines indicate the positions of 8 r ,m- The stress ratio is linear with tan 6, but flow will not occur unless the 
stress ratio remains above tan 6 r in the case of a system already flowing, or exceeds tan 9 m for initiation of flow. 



V. CONCLUSIONS 

We have performed a systematic simulation study of inclined plane, gravity-driven flow of granular materials, 
elucidating the reasons behind the discrepancies that exist between several different experimental procedures. These 
arise as a consequence of how far the system is from jamming, or how close the inclination angle of the pile is from 
the angle of repose. 

We have accurately determined the flow phase diagram for a system with a particular interaction parameter set, 
indicating the height dependence of the angle of repose and the maximum angle of stability. We have also determined 
the boundary between steady state and unstable flow and find that the maximum angle for steady state flow also 
exhibits height dependence. This suggests the existence of an intrinsic length scale in granular flows and forms part 
of our work in progress. 

For thick enough piles flowing on a rough inclined plane, say h > 20d, the velocity profiles and rheology follow 
Bagnold scaling. As one reduces the height of the flowing pile, a continuous transition from Bagnold rheology to 
linear velocity profiles through to avalanche- like dynamics occurs, until finally, one reaches the angle of repose 9 r , and 
flow ceases. As a result, one can obtain various velocity profiles and flow behaviors not only through changing the 
height of the flowing pile, but also by varying the inclination angle of the chute, such that there is an overlap region 



where one can obtain similar flow properties through either procedure. These features can be used in surface flow 
treatmentsilS to better predict the evolution of an avalanche surface. 

Our investigation on the initiation and cessation of flow shows that the initiation process starts from failure at the 
base that results in a "heat" front that propagates up the pile. Upon cessation of flow, an effective freezing front 
propagates up the pile with particles becoming trapped in the frozen layer until the whole system comes to rest. In 
some cases, flow initiation or cessation only occurred after many time units (approximately equivalent to 20 seconds 
for a 1mm sized particle). A consequence of these findings is that experimental reports of various flow behaviour 
in flowing piles may also be due to the time window of observation. In general, we have captured all the major 
observations from experiment, including quantitative agreement of the scaling parameter (3, Eq. ^, from Pouliquen's 
empirical theoryll. 
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